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Abstract 

The magnetic behavior of the low-dimensional phosphates (Sr,Ba)2Cu(P04)2 and BaCuP207 was 
investigated by means of magnetic susceptibility and ^^P nuclear magnetic resonance (NMR) mea- 
surements. We present here the NMR shift, the spin-lattice (1/Ti) and spin-spin (I/T2) relaxation- 
rate data over a wide temperature range 0.02 K < T < 300 K. The temperature dependence of 
the NMR shift K{T) is well described by the S = 1/2 Heisenberg antiferromagnetic chain model 
(D.C. Johnston et al, Phys. Rev. B 61, 9558 (2000)) with an intrachain exchange of J/ks — 165 
K, 151 K, and 108 K in Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuP207, respectively. Deviations 
from Johnston's expression are seen for all these compounds in the T-dependence of K{T) at low 
temperatures. '^^P is located symmetrically between the Cu ions and fluctuations of the staggered 
susceptibility at q = ^ should be filtered out due to vanishing of the geometrical form factor. 
However, the qualitative temperature dependence of our ^^P NMR 1/Ti indicates that relaxation 
due to fluctuations around Q = (but 7^ ^) have the same T-dependence as those at g = ^ and 
apparently dominate. Our measurements suggest the presence of magnetic ordering at 0.8 K in 
BaCuP207 {J/ks ~ 108 K) and a clear indication of a phase transition (divergence) in l/Ti(T), 
1/T2(T), and a change of the line shape is observed. This enables us to investigate the ID behavior 
over a wide temperature range. We find that 1/Ti is nearly T- independent at low-temperatures 
(1 K < T < 10 K), which is theoretically expected for ID chains when relaxation is dominated by 
fluctuations of the staggered susceptibility. At high temperatures, 1/Ti varies nearly linearly with 
temperature, which accounts for contribution of the uniform susceptibility. 

PACS numbers: 75.10.Pq, 75.40.Cx, 76.60.-k, 76.60.Cq 
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I. INTRODUCTION 



There is presently a lot of interest in the magnetic properties of one-dimensional (ID) 
Heisenberg antiferromagnetic (HAF) spin systems. This is due to the rich physics that 
they exhibit, in addition to the fact that such systems are tractable from a computa- 
tional/theoretical standpoint. In particular, qualitative differences are expected between 
integer-spin and half-integer spin HAF chains. While the integer-spin chains are gapped,— 
the half-integer spin chains are said to have quasi long-range order (LRO) due to the gradual 
spatial decay (power-law) of the spin-spin correlation function.i^ 

The magnetic Hamiltonian describing a spin-half Heisenberg chain can be written as 
H = — J Y^. Si-Si+i where J is the intrachain coupling constant between the nearest-neighbor 
spins. The temperature dependence of the magnetic susceptibility x(T) for the S = 1/2 
HAF chain was numerically calculated by Bonner and Fisher— and since then, the Bonner- 
Fisher expression has been used by experimentalists to determine the value of the exchange 
coupling (J) from the temperature dependence of the bulk susceptibility. A more accurate 
and analytical evaluation of the susceptibility of S* = 1/2 HAF chain was done by Eggert, 
Affleck, and Takahashi^ which is valid at low-temperatures. Via numerical simulations, an 
expression for x (T) accurate for both low and high temperatures (5 x 10^^^ < < 5, 
with ks the Boltzmann constant) was given by Johnston et al.^ . Dynamical properties of 
S = 1/2 chains have also been theoretically investigated. In particular, work has been 
focussed on the properties measured by nuclear magnetic resonance (NMR) techniques. 
Sachdev^ determined the temperature dependence of the NMR spin-lattice (l/Ti) and the 
Gaussian spin-spin (I/T2G) relaxation rates for half-integer spin chains for <^ 1 as: l/Ti 
= constant and I/T2G oc 1/\/T. Quantum Monte Carlo calculations by Sandvil# support 
these results over an appropriate temperature range. These results are at variance from 
those for classical spin chains {S = 00), where theory predicts- 1/Ti and I/T2G oc T~^/^. 

While the ID compounds CuCl2-2NC5H5and KCuFs have been experimentally investi- 
gated previously ,^°i^^ the onset of LRO (at the ordering temperature T/v) due to inter-chain 
interactions prevents the study of true ID properties down to low temperatures. For the 
abovementioned ID compounds, the ratios ^slk ^ 0.084 and 0.1195, respectively, have been 
determined. Dynamic and static properties of Sr2Cu03 ID chain have been extensively 
studied.iSi^'i^ii^ Due to the large value of J/ke (^ 2200 K) and the very weak inter-chain 
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couplings, this compound orders only below 5 K and hence ID properties could be studied 
in a large range of temperature. It is clearly useful to examine new compounds which might 
exhibit ID behavior in a large temperature range, thereby allowing for a comparison with 
theoretical models and improving our understanding of such systems. 

Sr2Cu(P04)2^^ and Ba2Cu(P04)2i^ are two isostructural compounds having a monoclinic 
unit cell with space group C2/m- The reported lattice constants are a = 11.515 A, 6 = 
5.075 A, c = 6.574 A and a = 12.160 A, 6 = 5.133 A, c = 6.885 A for Sr2Cu(P04)2 and 
Ba2Cu(P04)2, respectively. BaCuP207, which differs slightly in structure compared to the 
other two, crystallizes in a triclinic unit cell with space group PI and lattice constants 
a = 7.353 A, 6 = 7.578 A, c = 5.231 A.— In the former two compounds, each CUO4 
square plane shares it's edges with two similar kind of PO4 groups. The edge sharing takes 
place in one direction forming an isolated [Cu(P04)2]oo chain along the crystallographic b- 
direction. A likely interaction path of Cu^+ ions is sketched in Fig. 1(a). As opposed to 
Sr2Cu(P04)2 and Ba2Cu(P04)2, BaCuP207 contains two inequivalent P atoms, where each 
CUO4 plaquette shares it's edges with two different PO4 groups forming chains as shown in 
Fig. 1(b). Unlike the isolated chains of Sr2Cu(P04)2 and Ba2Cu(P04)2, there appear to be 
pairs of chains in BaCuP207. Detailed magnetic properties of these compounds have not 
been reported yet. Only Etheredge and Hwi>i^ have published the bulk susceptibility as a 
function of temperature for Ba2Cu(P04)2. However, the authors failed to comment on the 
broad maximum at 80 K presumably because it was suppressed by a high Curie contribution 
present in their sample. 

In this paper we present in detail the magnetic properties of the ID copper phosphates, 
Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuP207 using ^^^P NMR as a local probe. NMR is 
regarded as a valuable tool for the study of microscopic properties of ID chains, especially 
through the studies of the NMR shift (K), the spin-spin relaxation rate (I/T2), and the spin- 
lattice relaxation rate (1/Ti). We report on measurements of the bulk susceptibility x{T) 
for 1.8 K < T < 400 K and K{T), l/Ti(T), and 1/T2(T) of ^^P NMR in a large temperature 
range (0.02 K < T < 300 K). The experimental details concerning sample preparation and 
various measurements are given in the next section. Section III contains our experimental 
results and a detailed discussion about the results is presented in section IV. Our work on 
these compounds strongly suggests that they are some of the best examples of 5* = 1/2 
ID HAF systems. In the course of our work, magnetic susceptibility and heat capacity of 
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Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuP207 were reported by Belik et al^^^. They found 
the exchange constant {J/ks) to be 144 K for Sr2Cu(P04)2, 132 K for Ba2Cu(P04)2, and 
103.6 K for BaCuP207. Presence of any LRO was not seen from specific heat measurement 
down to 0.45 K for Sr2Cu(P04)2 and Ba2Cu(P04)2, whereas BaCuP207 showed ordering at 
0.81 K. 

II. EXPERIMENTAL DETAILS 

Polycrystalhne samples of Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuP207 were prepared 
by sohd state reaction techniques using BaCOs (99.9% pure), SrCOs (99.999% pure), CuO 
(99.99% pure) and (NH4)2HP04 (99.9% pure) as starting materials. The stoichiometric mix- 
tures were fired at 800 °C (Sr2Cu(P04)2), 700 °C (Ba2Cu(P04)2) and 650 °C (BaCuP207) 
for 120 hours each, in air, with several intermediate grindings and pelletization. Finally 
some amounts of each of the samples were annealed at 400 °C under a reducing atmosphere 
(5% H2 in Ar) in an attempt to reduce the Curie contribution in the bulk susceptibility. 
Nearly single phases were confirmed from x-ray diffraction, which was performed with a 
Philips Xpert-Pro powder diffractometer. A Cu target was used in the diffractometer with 
Xav = 1.54182 A. The x-ray diffraction patterns are shown in Fig. 2. Only one compound 
(Sr2Cu(P04)2) exhibits minor impurity peaks which are marked by asterisks. Lattice pa- 
rameters were calculated using a least-square fit procedure. The obtained lattice constants 
are (a = 11.496(5) A, b = 5.069(2) A, c = 6.566(3) A), (a = 12.138(2) A, b = 5.123(1) A, 
c = 6.868(1) A) and (a = 7.338(2) A, 6 = 7.561(2) A, c = 5.217(1) A) for Sr2Cu(P04)2, 
Ba2Cu(P04)2, and BaCuP207, respectively. These are in agreement with previously re- 
ported values. 

Magnetization (M) data were measured as a function of temperature T (1.8 K < T < 
400 K) and applied field H {0 < H < 70 kG) using a SQUID magnetometer. 

The NMR measurements were carried out using pulsed NMR techniques on ^^P nuclei 
(nuclear spin 1 = 1/2 and gyromagnetic ratio 7/27r = 17.237 MHz/Tesla) in a large tem- 
perature range (0.02 K < T < 300 K). We have done the measurements at two different 
applied fields of about 55 kG and 4 kG which correspond to radio frequencies (rf) of about 
95 MHz and 6.8 MHz, respectively. 

For 2 K < T < 300 K, NMR measurements were done in a 55 kG applied field with a 
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^He cryostat (Oxford Instruments). Spectra were obtained by Fourier transform (FT) of 
the NMR echo signals using a 7r/2 pulse of width of about 4 /is. The NMR shift K{T) = 
[z^ (T) — Uref] /j^ref was determined by measuring the resonance frequency of the sample 
(z/(T)) with respect to a standard H3PO4 solution (resonance frequency z^re/)- The spin- 
lattice relaxation rate (1/Ti) was determined by the inversion-recovery method. Spin-spin 
relaxation rate {I/T2) was obtained by measuring the decay of the echo integral with variable 
spacing between the 7r/2 and the tt pulse. 

In the 0.02 K < T < 10 K range, NMR measurements were performed using a ^He/^He 
dilution refrigerator (Oxford Instruments) with the resonant circuit inside the mixing cham- 
ber. Spectra were obtained by field sweeps at a constant radio frequency (z^r/) of 95 MHz. 
1/Ti was measured down to 0.02 K following the same procedure as described above using 
7r/2 pulse of width 15 /is. Lower rf power (and consequently longer pulse widths) were used 
to avoid rf heating of the sample. Measurements were also done in a low-field of about 4 kG 
{vrf ~ 6.8 MHz) where the NMR line was narrow and inversion of the nuclear magnetiza- 
tion by a TT pulse of width 30 /is was assured. The data from low-field measurements almost 
reproduce the high-field data. 

III. RESULTS 

A. Bulk Susceptibility 

Magnetic susceptibilities x(r) (= M/H) for all the three compounds were measured as 
a function of temperature in an apphed field of 5 kG (Fig. 3). The amount of ferromagnetic 
impurities present in our samples were estimated from the intercept of M vs. H isotherms 
at various temperatures and were found to be 19 ppm, 12 ppm, and 30 ppm of ferromagnetic 
Fe'^^ ions for Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuP207, respectively. The data in Fig. 3 
have been corrected for these ferromagnetic impurities. As shown in the figure, all the sam- 
ples exhibit a broad maximum, indicative of low-dimensional magnetic interactions. With 
decrease in temperature, susceptibility increases in a Curie- Weiss manner. This possibly 
comes from chain ends, natural defects, excess oxygen and extrinsic paramagnetic impuri- 
ties present in the samples. No obvious features associated with LRO are seen for 1.8 K < 
T < 400 K for any of the samples. A substantial reduction of Curie terms was achieved by 
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annealing the samples at 400 ''C in an atmosphere of 5% H2 in Ar. Similar experiments in 
Sr2Cu03 and Y2BaNi05 lead to reduced Curie terms .^^i^^ In the insets of Figs. 3(a) and 
3(b), data for as-prepared and reduced Sr2Cu(P04)2 and Ba2Cu(P04)2, respectively are 
shown. Since the Curie contribution in the case of as-prepared BaCuP2 07 is not large, we 
did not treat this sample in a reducing atmosphere. 

In order to fit the bulk susceptibility data, we assume that the susceptibility consists of 
three terms: 

C 

X = Xo + 7777-^ + Xspm(T) (1) 



where Xspin{T) is the uniform spin susceptibility for aS' = l/2 1D HAF system given in 
Ref. 6 (expression corresponding to "fit2"). This expression (containing the the Lande g- 
factor and J as fitting parameters) is not reproduced here since it is somewhat unwieldy. The 
first term Xo is temperature independent and consists of diamagnetism of the core electron 
shells (xcore) and Van-Vleck paramagnetism (xw) of the open shells of the Cu^"*" ions present 
in the sample. The second term is the Curie- Weiss contribution due to paramagnetic 
species in the sample. 

The average Lande (/-factors determined from an analysis of the powder spectra from 
electron paramagnetic resonance (EPR) experiments on our samples were found to be 2.15, 
2.15, and 2.2 for Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuP207, respectively. Our experimen- 
tal x(^) data were fitted using the above (/-values (the solid lines are the best fits in Fig. 3) 
and the extracted parameters are listed in Table I. 



Table I. 



Sample 


Xo 


C ^ 




10"^ cm^/mole 10" 


cm^K/mole K K 


Sr2Cu(P04)2 


0.005 


10.6 1.1 152 


Ba2Cu(P04)2 


-0.15 


6.8 0.5 146 


BaCuPaOr 


-0.07 


1.6 0.4 108 



Adding the core diamagnetic susceptibility for the individual ionsi2^, the total Xcore was 
calculated to be — 1.39x 10""^ cm^/mole, — 1.73x 10""^ cm^/mole, and — 1.29x 10^'^ cm^/mole 
for Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuP207, respectively. The Van-Vleck paramagnetic 
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susceptibility for our samples estimated by subtracting Xcore from xo gives Xw = 14.4 x 10~^ 
cm^/mole, 2.3 x 10~^ cm^/mole, and 5.9 x 10~^ cm^/mole for Sr2Cu(P04)2, Ba2Cu(P04)2, 
and BaCuP207, respectively. These values are comparable to that found in Sr2Cu03^. The 
Curie contributions present in the samples correspond to a defect spin concentration of 3 
%, 1.8 % and 0.4 % for Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuP207, respectively assuming 
defect spin S = 1/2. 

B. 3ip NMR 

1. NMR shift 

NMR has an important advantage over bulk susceptibility for the determination of mag- 
netic parameters. While the presence of a Curie-like contribution restricts the accurate 
determination of Xspin{T) and hence J from x{T)i in NMR this paramagnetism broadens 
the NMR line but does not contribute to the NMR shift K. Therefore, it is more reliable to 
extract the Xspin{T) and J from the temperature dependence of the NMR shift rather than 
from the bulk susceptibility. From Fig.l it appears that in all the compounds each '^^P is 
coupled to two Cu^^ ions via a super-transferred hyperfine interaction mediated by oxygen 
ions in its neighborhood. All the NMR data reported in this paper correspond to samples of 
Sr2Cu(P04)2 and Ba2Cu(P04)2 which were treated in a reducing atmosphere as described 
in the previous section while the data for BaCuP207 correspond to the as-prepared sample. 
We note here that we also did NMR measurements on the as-prepared Sr2Cu(P04)2 and 
Ba2Cu(P04)2 samples (for T > 10 K) and found no differences with respect to the samples 
which were treated in a reducing atmosphere. 

NMR shift data as a function of temperature are shown in Fig. 4. The samples ex- 
hibit broad maxima at different temperatures: ~ 100 K for Sr2Cu(P04)2, ~ 90 K for 
Ba2Cu(P04)2, and 70 K for BaCuP207, indicative of short-range ordering. Towards 
lower temperatures T < 20 K, the NMR shift K{T) shows a plateau as is demonstrated by 
the semilogarithmic plot in the inset of Fig. 4. In the sub-Kelvin region, the NMR shift 
K{T) of all our samples decreases steeply. The fall-off appears below ksT / J ~ 0.003 for 
Sr2Cu(P04)2, keT/J ~ 0.0033 for Ba2Cu(P04)2, and keT/J ~ 0.01 for BaCuP207. 

The conventional scheme of analysis is to first determine the spin susceptibility (as done 
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in the previous section) and then plot K vs XSpin with T as an imphcit parameter. The 
exchange couphng J is obtained from the susceptibihty analysis while the slope of the K 
vs XSpin plot yields the total hyperfine coupling A between the ^^P nucleus and the two 
nearest-neighbor Cu^+ ions. As an example, such a. K — Xspin plot is shown in Fig. 5 for 
BaCuPaOy. 

Since an algebraic expression for the temperature dependence of the spin susceptibil- 
ity (and therefore the spin-shift) is known in this case, we prefer to determine J and A 
simultaneously by fitting the temperature dependence of K to the following equation, 

where Kq is the chemical shift and Na is the Avogadro number. While fitting, g was 
kept fixed to the value obtained from EPR analysis and Kq, A, and J were free parameters. 
The parameters J and A determined in this manner arc considered more reliable, since the 
only temperature dependent term in the NMR shift is due to spin-susceptibility, while bulk 
susceptibility analysis is somewhat hampered by low-temperature Curie terms. As shown 
in Fig. 4, the shift data fit nicely to Eq. 2 in the temperature range 10 K < T < 300 K 
yielding the parameters given in Table II. 



Tabic II. 



Sample 


Ko 


A J/ks 




ppm 


Oc//is K 


Sr2Cu(P04)2 


47 


2609 ± 100 165 


Ba2Cu(P04)2 


40 


3364 ± 130 151 


BaCuPaOv 


73 


2182 ±20 108 



2. Spectra 

For all the three compounds the ^^P NMR spectra consist of a single spectral hne as 
is expected for / — 1/2 nuclei (Fig. 6). As shown in the crystal structures in Fig. 1, 

Sr2Cu(P04)2 and Ba2Cu(P04)2 have a unique ^^P site, whereas in BaCuP207 there are 
two incquivalcnt ^^P sites. However, a single resonance line even for BaCuP207 implies 
that both the ^^P sites in this compound are nearly identical. Since our measurements are 
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on randomly oriented polycrystalline samples, asymmetric shape of the spectra corresponds 
to a powder pattern due to an asymmetric hyperfine coupling constant and an anisotropic 
susceptibility. The linewidth was found to be field and temperature dependent as is shown 
in the insets of Figs. 6(a) and (b) for (Sr/Ba)2Cu(P04)2 and in Figs. 6(c) and (d) for 
BaCuP207, respectively. While we did not do a detailed analysis of the linewidth, its T- and 
iJ-dependence is likely due to macroscopic field inhomogeneities due to the demagnetization 
effects of a powder sample^^ and paramagnetic impurities. 

As seen from Fig. 6 (c) the NMR spectra of BaCuP207 broaden abruptly below about 0.8 
K. We then measured the spectral lineshape in a low field {H ^ 4 kG) below 0.8 K in order 
to check whether any features could be resolved. Fig. 6 (d) shows the appearance of two 
shoulders on either side of the central line below 0.8 K. This is most likely an indication of the 
appearance of LRO. The positions of the shoulders stay unchanged with temperature while 
their relative intensity increases with decreasing temperature. A more detailed discussion is 
carried out in section IV. 



3. Spin-lattice relaxation rate l/Ti 

Temperature dependencies of ^^P 1/Ti are presented in Fig. 7. For the 1/Ti experiment, 
the central positions of corresponding spectra at high (55 kG) and low (4 kG) external fields 
have been irradiated. Inset of Fig. 7 (a) shows the typical magnetization recovery at i7 ~ 
55 kG and at two different temperatures. For a spin-1/2 nucleus the recovery is expected to 
follow a single exponential behavior. In Sr2Cu(P04)2 and Ba2Cu(P04)2 (for if ^ 55 kG), 
the recovery of the nuclear magnetization after an inverting pulse was single exponential 
down to 2 K while for T < 2 K, it fitted well to the double exponential, 

' ^ -^lexp -— +A2exp -— +C (3) 



2 V M(oo) ; ^ " V T,lJ V T,s, 

where 1/Til corresponds to the slower rate and I/T15 is the faster component. M{t) is the 
nuclear magnetization a time t after an inverting pulse. Since the deviation from single 
exponential behavior could be due to the large linewidth and our consequent inability to 
saturate the NMR line, we also performed experiments at a lower field (~ 4 kG), where the 
line is about three times narrower (see insets of Figs. 6 (a) or (b)). However, even in low- field 
(where the rf field Hi was sufficient to ensure complete inversion) the nuclear magnetization 
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recovery is not single exponential implying that this is an intrinsic effect. With increasing 
temperature, the ratio ^ decreases and the recovery becomes single exponential for T > 
2 K. It appears that the longer Ti component comes from the chain itself while the faster 
component is associated with ^^P nuclei near chain ends. Clearly, at lower temperatures, 
the chain-end- induced magnetization extends to large distances from chain ends (thereby 
affecting more ^^P nuclei) and consequently the weight associated with the faster relaxation is 
more at lower temperatures. From the experiment, it was found that low-field measurements 
reproduce almost the same Ti as for high-field.— Although we have measured Ti down to 0.02 
K, since the magnetization recovery is not single-exponential, reliable relaxation rates l/Ti 
couldn't be obtained below 0.1 K (where the faster component l/Ti^ has a large weight). 
Fig. 7 (a) and (b) display data down to 0.1 K, where it is seen that 1/Ti for Sr2Cu(P04)2 
and Ba2Cu(P04)2 do not show any anomaly. Even at lower temperatures, there was no 
indication of a divergence of the relaxation rate, indicating the absence of any magnetic 
ordering. For 1 K < T < 10 K, 1/Ti remains constant with temperature and below 0.5 K a 
slight increase was observed for both Sr2Cu(P04)2 and Ba2Cu(P04)2. At high temperatures 
(T > 30 K), 1/Ti varies nearly linearly with temperature. 

In BaCuP207, 1/Ti at ~ 55 kG was measured down to 3 K. Once again, the large 
line width prevented us from saturating the nuclear magnetization below 3 K. Low-field 
measurements give perfect single exponential recovery down to 2 K and below 2 K, it was 
fitted well to double exponential. From Fig. 7(c), it is clear that the l/Ti(T) diverges at T 
~ 0.8 K, indicating an approach to magnetic ordering. For 1 K < T < 10 K, 1/Ti remains 
constant and for T > 15 K, it varies linearly with temperature. Slight change in magnitude 
in low-field data compared to high-field data may be due to spin-diffusion.i^^ 

4- Spin-spin relaxation rate I/T2G 

Spin-spin relaxation was measured at if ^ 4 kG, where the line is sufficiently narrow. 
The spin-spin relaxation was found to have a Gaussian behavior and the rate (I/T2G) was 
obtained by monitoring the decay of the transverse magnetization after a | — t — vr pulse 
sequence, as a function of the pulse separation time t, and fitting to the following equation, 

M {2t) = Mo exp[-2 (^^^ ]+C (4) 
10 



As shown in the inset of Fig. 8, the spin-spin relaxation rate I/T2G for all the samples 
is nearly temperature independent. BaCuP207 is the only compound which exhibits a 
significantly enhanced spin-spin relaxation rate I/T2G at the lowest temperature compared 
to elevated temperatures. This increase of I/T2G in BaCuP207 is most likely related to 
LRO. 

IV. DISCUSSION 
A. NMR shift 

The general variation of the shift with temperature follows the expected behavior of 
an S = 1/2 HAF chain, as seen in the results section. A steep decrease in K{T) was 
observed below fcijT/J ~ 0.003 for Sr2Cu(P04)2, keT/J ~ 0.0033 for Ba2Cu(P04)2 and 
ksT/ J ~ 0.01 for BaCuP207.^^ This decrease of K is clearly much more than the logarithmic 
decrease with an infinite slope at zero temperature, expected from theory (see solid line in 
Fig. 4). In Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuP207, from Fig. 4, the extrapolated shifts 
at zero temperature were found to be 400 ppm, 590 ppm, and 580 ppm, respectively. The 
theoretically expected values, derived from the T = susceptibility^^ and using the relevant 
A and Kq, are 544 ppm, 740 ppm, and 738 ppm, respectively for Sr2Cu(P04)2, Ba2Cu(P04)2, 
and BaCuP207. Among the various causes for this deviation, one might be the onset of spin- 
Peierls order. In such a case, the spin susceptibility (and therefore the spin shift) should 
vanish at T = 0. However, our extrapolated T = shifts are much more than the chemical 
shifts Ko and there is no exponential decrease of l/Ti(T) toward low temperatures. Another 
possibility is the onset of 3D LRO. In this case, a divergence should have been seen in the 
temperature dependencies of the spin-lattice relaxation rate 1/Ti as well as in the spin-spin 
relaxation rate 1/T2G. While this is the case for BaCuP207, only a small increase of the 
relaxation rates is observed for Sr2Cu(P04)2 and Ba2Cu(P04)2. A clear effect is observed 
in the temperature dependencies of K(T), 1/Ti, I/T2, and lineshape for BaCuP207 at 0.8 
K. This establishes the onset of LRO at 0.8 K in BaCuP207. However, in the case of 
Sr2Cu(P04)2 and Ba2Cu(P04)2, while a clear anomaly is seen in K{T) at low-temperature, 
only a weak anomaly is seen in l/Ti(T) and no significant changes were observed either in 
the low-temperature spectra or in 1/T2(T). In summary, the presence or absence of LRO at 
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low-temperatures in Sr2Cu(P04)2 and Ba2Cu(P04)2, can't be unambiguously established. 



B. Effect of chain ends on the NMR spectrum 

From field-theory and Monte Carlo calculations, Eggert and Affleck— found that in case of 
half-integer spin chains, the local susceptibility near an open end of a finite chain has a large 
alternating component. This component appears in the form of staggered magnetization 
near chain ends, under the influence of an uniform field. This staggered moment has a 
maximum at a finite distance from the end and increases as 1/T with decreasing temperature. 
Analytical expression for the spin susceptibility x (0 site I were obtained which consist 
of the uniform and alternating {Xait) parts^^ x (0 = Xu (0 + ("1)' Xait (0 ^ where the 
uniform part is nearly constant and the alternating part is given by 

= V^(^;/7rT)sinh {27(Tl/v) 



where v = 7rJ/2 is the spin- wave velocity. The NMR spectrum represents the dis- 
tribution function of NMR shift, which is equivalent to Xait, and has the form g (x) = 



fx — (—1) Xait (0 5 where / takes the form of a Lorentzian. This expression is valid 
for in-chain Cu^^ site. In our compounds, ^^P is the probe nucleus, which is sensitive to two 
nearest neighbor Cu^+ ions belonging to one chain. ^^P NMR line shape is then given by 
9{x) =Y.if ^ - {(-!)' Xait (0 + (-1)'^^ Xait + 1)} • Fig. 9 shows the simulated spectra 
for both Cu and P sites choosing / to be a Lorentzian with width 0.05. 

For ksT = J/33, at Cu site, Xait has a maximum at I = 0.48J/T, which results in features 
in the spectra on either side of the central line. With increasing temperature these features 
on either side of the central peak disappear and VT{-^^) (where AH is the width at the 
background and Hq is the field at the central peak) remain constant with temperature. This 
has been seen by Takigawap^ in ^^Cu NMR spectra of Sr2Cu03. However, the static effects 
of the staggered magnetization at the ^^P site in our compounds are expected to be much 
weaker (Fig. 9) due to the near cancellation of the magnetization from the neighboring Cu 
sites. 

Our low-temperature spectra for Sr2Cu(P04)2 and Ba2Cu(P04)2 show a single spectral 
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line without any shoulders on either sides of the central peak. This is in agreement with 
our expectations that we are not sensitive to chain-end effects in ^^P NMR spectra of these 
compounds. 

In low-field NMR on BaCuP207, a sudden increase of linewidth was observed below 
0.8 K along with appearance of two shoulder-like features located symmetrically on either 
side of the central peak (Fig. 6(d)). If these shoulder-like features come from the staggered 
magnetization of chain ends then, as discussed above, VT{j^) should be temperature in- 
dependent with the shoulders moving outwards at lower temperatures. But in BaCuP207 
the shoulder positions are temperature independent, suggesting that those features are not 
chain-end effects. Further, in case of a structural phase transition, symmetrically located 
features which become more intense as temperature is lowered are not expected. Our data 
are consistent with an alternating magnetization in the Cu chain. Further, the magnitude 
of this magnetization must have a maximum away from the chain end (so as to produce the 
symmetric features) and the location (site Imax) of the maximum must not change with T 
(as the shoulder positions are temperature independent). As shown before, the chain-end 
induced effects are expected to be negligible at ^^P site. The large effect seen here must 
stem from a magnetic transition possibly with exotic spin-order. 

C. Wave vector q- and temperature T- dependence of 1/Ti 

In order to study the microscopic behavior of ID HAF systems, it is useful to measure the 
temperature dependence of spin-lattice relaxation rate which yields information on the imag- 
inary part of the dynamic susceptibility x (Qj ^)- The spin-lattice relaxation rate, in general, 
is affected by both uniform (g = 0) and staggered spin fluctuations [q = ±^)- The uniform 
component leads to 1/Ti oc T, while the staggered component gives 1/Ti = constant.^ At 
the ^^P sites, g-dependence of 1/Ti can be expressed in terms of form factors as. 



We are probing on the ^^P nucleus, which is linked to the Cu spins via 0(1) and 0(2) 
as shown in Fig. 1. Since ^^P is symmetrically located between the Cu ions, the antiferro- 
magnetic fluctuations are filtered at the ^^P site, provided the two hyperfine couplings are 
equal. However, contributions just slightly differing from q = n/a are apparently strong 




(6) 



g 
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enough and/or the two hyperfine couphngs are unequal and result in a qualitative behavior 
expected when relaxation is dominated by fluctuations of the staggered susceptibility. 

If one were to ignore the geometrical form factor completely, the relaxation rate due to 
staggered fluctuations can be calculated following the prescription of Barzykin^^. He ob- 
tained the normalized dimensionless NMR spin-lattice relaxation rate at low-temperature 
(l/Ti) = .^"L 0.3, where Ath is A{2h'j/27r). Assuming the fluctuations to be cor- 



related, 1/Ti can be written as 1/Ti = tjta2- Using this expression, (1/Ti) at the '^^P 



th 

rio ixm+'^on qc 1 / # - — 

hJ/A^ 

site was calculated to be about 44 sec,^^ 80 sec,~^ and 47 sec,^^ whereas our experimental 
values are 15 sec,~^ 20 sec,~^ and 25 sec,~^ for Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuP207 
respectively in the 1 K < T < 10 K range. The experimental values are clearly smaller than 
the theoretical ones due to the geometrical form factor. Further, a logarithmic increase of 
1/Ti is expected at low temperatures followingi^ 




(7) 



where D = 1/(27?)^/^, A is the cutoff parameter given by 2v^e'^+V and C {c^ 0.5772157) 
is Euler's constant. In Fig. 10, we have plotted {{i/Ti)/{l/Ti)T=ioK)norm calculated from 
equation 7 and the experimental results of the spin-lattice relaxation rate 1/Ti for our 
compounds normalised by their values at 10 K, as a function of temperature. The quali- 
tative temperature dependence of the experimental spin-lattice relaxation rate 1/Ti agrees 
reasonably well with theoretical calculations in the constant region but below 0.5 K, the 
experimental increase in Sr2Cu(P04)2 and Ba2Cu(P04)2 is somewhat more than the loga- 
rithmic increase expected theoretically. However, complementary measurements are needed 
to fully understand this issue. 



D. Spin-spin relaxation rate I/T2G 

Following the treatment of Sachdevi and Barzykin,— spin-spin relaxation is expected 
to follow a Gaussian behavior with a temperature dependence given by I/T2G' oc \/\/T. 
Further, the spinon mediated spin-spin relaxation rate can be calculated^ as follows. The 
normalized spin-spin relaxation rate is given by (^f^j = (^^) a^t ■ Dividing 
ij^l by (^) and equating it to 1.8 (the value obtained from Ref. 16) one finds 

V ^ / norm V 2G / norm 
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/ \ 1/2 / \ 1/2 / \l/2 

that I/T2G = ^ (^) , H (^) , and § , respectively for Sr2Cu(P04)2, 

Ba2Cu(P04)2, and BaCuP207. This leads to T2G values 3.2 ms, 1.8 ms, and 3.6 ms in 
contrast to our experimental values 269 us, 255 us, and 207 us respectively for Sr2Cu(P04)2, 
Ba2Cu(P04)2, and BaCuP207 at 1 K. These are almost three orders of magnitude smaller 
than the theoretically calculated values. Further, our experimental spin-spin relaxation 
rates are temperature independent. Clearly, in the present case, spinon mediated coupling 
does not contribute to spin-spin relaxation. On the other hand, an estimate of the nuclear 
dipole-dipole mediated relaxation = where r is the dipole-dipole distance^ would 
seem to explain the observed relaxation rates. This must be primarily because of the small 
exchange coupling in contrast to Sr2Cu03 where J is an order of magnitude larger and 
spinon mediated Gaussian spin-spin relaxation rate has been observed. 

V. CONCLUSION 

Our NMR and susceptibility measurements on Sr2Cu(P04)2, Ba2Cu(P04)2, and 
BaCuP207 show good agreement with the theory of ID S* = 1/2 Heisenberg antiferromag- 
netic chains. NMR shift Kas a function of temperature fitted well to the recent theoretical 
calculation by Johnston and the exchange interaction J/ks is estimated to be 165 K, 151 
K, and 108 K for Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuPaOy, respectively. We observed 
a steep decrease of the NMR shift K{T) below T ~ 0.003 J/ A;^, 0.0033 J/Zcs, and OMJ/ke 
for Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuP207, respectively. Low-field ^^P NMR spectra 
of BaCuP207 shows sudden appearance of broad humps on either side of the central peak 
for T < 0.8 K, indicating the onset of LRO. The spin-lattice relaxation rate 1/Ti was mea- 
sured in a temperature range 0.02 K < T < 300 K. No clear indication of any kind of 
magnetic ordering was seen in Sr2Cu(P04)2 and Ba2Cu(P04)2 from the 1/Ti data, whereas 
a clear indication of magnetic ordering was observed at T f=:i 0.8 K (/csT/J 0.0074) for 
BaCuP207. At low-temperature, 1/Ti follows a nearly logarithmic increase for T < 0.5 
K for Sr2Cu(P04)2 and Ba2Cu(P04)2, which is expected for a ID S* = 1/2 HAF system. 
Though the transverse decays follow Gaussian behavior for all our samples, they result from 
dipole-dipole interaction rather than a spinon mediated interaction. Our experimental evi- 
dences on these compounds strongly reflect their low-dimensional nature making them one 
of the best ID 5" = 1/2 HAF systems that have been looked at so far. 
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Figure Captions 

FIG. 1 Schematic diagram of [Cu(P04)2]j^ linear chains (a) propagating along 6-direction 
for (Ba/Sr)2Cu(P04)2 and (b) along the c-direction for BaCuP207. The arrows show the 
direction of chains while the possible interaction paths are shown by solid lines between 
atoms. 

FIG. 2 X-ray diffraction pattern (Intensity vs. 26) for (a) Sr2Cu(P04)2, (b) Ba2Cu(P04)2, 
and (c) BaCuP207. Miller indices are not shown in order to maintain the clarity of the 
figures. Impurity peaks are marked by asterisks. 

FIG. 3 Magnetic susceptibility (M/H) vs. temperature T for (a) Sr2Cu(P04)2, (b) 
Ba2Cu(P04)2, and (c) BaCuP207 in an applied field of 5 kG. The solid lines are best fits of 
the data to Eq. 1. Insets of (a) and (b) display the comparison of the magnetic susceptibility 
for as-prepared and reduced samples at lower temperatures. Susceptibility variations as a 
function of ^ are shown to demonstrate the lowering of Curie term on reduction. 



17 



FIG. 4 shift K vs. temperature T for (a) Sr2Cu(P04)2, (b) Ba2Cu(P04)2, and (c) 
BaCuP207. The sohd hnes are fits of Eq. 2 in the temperature range, 10 K < T < 300 K 
and then extrapolated down to 0.01 K. Inset shows K vs. T on a logarithmic temperature 
scale for improved visualisation of the low-T data. 

FIG. 5 ^^P shift K vs. spin susceptibility Xspin for BaCuP207. The solid line shows the 
linear fit. 

FIG. 6 Low-field (// ~ 4 kG) ^^P NMR spectra at different temperatures T for (a) 
Sr2Cu(P04)2 and (b) Ba2Cu(P04)2. The insets of (a) and (b) contain the high-field and 
low-field spectra at 1 K. (c) and (d) show high-field and low-field spectra, respectively, for 
BaCuP207 at various temperatures below 1 K, showing the sudden change in line width 
there. 

FIG. 7 Spin-lattice relaxation rate l/Ti (both high and low fields) vs. temperature T for 
(a) Sr2Cu(P04)2, (b) Ba2Cu(P04)2, and (c) BaCuPaOy. In the inset of (a), the high-field 
magnetization recoveries are plotted as a function of pulse separation t and the solid line is a 
single-exponential fit for Sr2Cu(P04)2. The inset of (b) displays the relaxation rate data for 
Ba2Cu(P04)2 on a linear temperature scale. In the inset of (c), the low-field magnetization 
recoveries arc plotted as a function of pulse separation t at 1 K and 2 K. The 1 K data are 
fitted to double exponential (Eq. 3), while 2 K data arc fitted to single exponential. 

FIG. 8 Spin-echo decays are plotted as a function of t'^ at two different temperatures for 
(a) Sr2Cu(P04)2, (b) Ba2Cu(P04)2, and (c) BaCuPaOy. The sohd hnes show the fitting to 
a Gaussian function (Eq. 4). In the insets, I/T2G is plotted as a function of temperature T. 

FIG. 9 Distribution function g{x) which represents the NMR spectrum is plotted for 
(a) Cu site and (b) P sites. In the insets of (a) and (b), (— l)'xa/t(0 (~l)'xa«t(0 + 
(—1)'^^ Xait (l + 1) respectively are plotted as a function of the site index / from the chain 
end at two different temperatures. 

FIG. 10 Normahzed experimental data (open symbols) as well as theoretical curve (solid 
line) of {l/Ti)norm are plotted vs. ksT/J for Sr2Cu(P04)2, Ba2Cu(P04)2, and BaCuPaOy 
at low temperatures, 0.1 K < T < 10 K. The data and the theoretical curves have been 
scaled to 1 at 10 K. 

Table Captions 
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TABLE I. Values of the parameters (xo, C*, d-, and J//cb) obtained by fitting the bulk 
susceptibility to Eq. 1 to for each of the three samples. 

TABLE IL Values of the parameters {Kq, A, and J/Zcs) obtained by fitting the NMR 
shift data to Eq. 2 in the temperature range, 10 K < T < 300 K for each of the three 
samples. 
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